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Abstract 

This article presents recent work to constrain the physical and chemical 
properties in high-mass star formation based largely on interferometric 
high-spatial-resolution continuum and spectral line studies at (sub)mm 
wavelengths. After outlining the concepts, potential observational tests, 
a proposed evolutionary sequence and different possible definitions for 
massive protostars, four particular topics are highlighted: (a) What are 
the physical conditions at the onset of massive star formation? (b) 
What are the characteristics of potential massive accretion disks and 
what do they tell us about massive star formation in general? (c) How 
do massive clumps fragment, and what does it imply to high-mass star 
formation? (d) What do we learn from imaging spectral line surveys 
with respect to the chemistry itself as well as for utilizing molecules as 
tools for astrophysical investigations? 

1 Introduction 

Star formation is a key process in the universe, shaping the structure of en- 
tire galaxies and determining the route to planet formation. In particular, 
the formation of massive stars impacts the dynamical, thermal and chemical 
structure of the interstellar medium (ISM), and it is almost the only mode of 
star formation observable in extragalactic systems. During their early forma- 
tion, massive star-forming regions inject energy to the ISM via their outflows 
and jets, then during their main sequence evolution, the intense uv-radiation 
of the massive stars heats up their environment, and at the very end of their 
life, Supernovae stir up the ISM by their strong explosions. Furthermore, 
massive stars are the cradles of the heavy elements. Hence, life as we know 
it today would not exist if massive stars had not formed flrst the heavy el- 
•'-H . ements in their interiors via nucleosynthesis. In addition to this, almost all 
stars form in clusters, and within the clusters massive stars dominate the 
. overall luminosities. Isolated low-mass star formation is the exception, and 
■ " " ' isolated high-mass star formation likely does not exist 

Concepts: In spite of their importance, many physical processes during 
the formation of massive stars are not well understood. While there ex- 
ists a paradigm for low-mass s tar formation on which large parts of the sc i- 
entific community agree (e.g., Andre et al. 20Q(i McKee fc Qstriker 2007 ). 
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this is less the case for high-mass star formation (e.g., Beuther et al. 2007at 



Zinnecker &: Yorke.2007f ) . The conceptional problem is based on the fact that 



at least in spherical symmetry the radiation pressure of a centrally ignited 
star >8M0 would be large enough to stop any further gas infall and hence 
inhibit the formation of more massive objects. Over the last decade, two 
schools have been followed to solve this problem: (a) The turbulent accre- 
tion scenario, which is largely an enhancement of the low-mass star forma- 
tion scenario, forms massive stars in a turbulent core with high accretion 
rates and a geometry including accretion disks and molecular outflows (e.g., 
Yorke fc Sonnhaltei]l2Qoi iMcKee fc Ta^l20Q3t iKrumholz et al.ll2007bt iKetd 



2007l l. In contrast to that, (b) the competitive accretion scenario relies on 
the clustered mode of massive star formation. In this scenario, the accretion 
rates are determined by the whole cluster potential, and those sources sitting 
closest to the potential w ell will competitively accrete most of the mass (e.g.. 



Bonnell et al.ll20Q4l . 120071 1. 



Potential observational tests: How can we observationally discriminate 
between the different scenarios? For example, the turbulent accretion sce- 
nario predicts qualitatively similar outflow and disk properties as known for 
low-mass stars, however, with quantitatively enhanced parameters like the 
accretion rates, outflow energies, or disk sizes. In contrast to that, modeling 
of the proto-cluster evolution in the competitive accretion scenario indicates 
extremely dynamic movements of all cluster-members throughout their whole 
evolution. It is unlikely that in such a dynamic environment collimated out- 
flows or large disks could survive at all. Another difference between both 
scenarios is based on their early fragmentation predictions. While the turbu- 
lent accretion scenario d raws their initial gas clump s from (gravo-)turbulently 
fragmented clouds (e.g., Padoan fc Nordlundl 2002 ) and these gas clumps do 
not fragment afterwards much anymore, the competitive accretion scenario 
predicts that the initial gas clumps fragment down to many clumps all of 
the order a Jeans-mass (~O.5M0). Hence, while in the former scenario, the 
Initial Mass Function (IMF) is determined during the early cloud fragmenta- 
tion processes, the latter models predict that the IMF only develops during 
the ongoing cluster formation. Therefore, studying the initial fragmentation 
and the early core mass functions can give insights in the actual massive star 
formation processes. 

Evolutionary sequence: Independent of the formation scenario, there has to 
be an evolutionary sequence in which the processes tak e place. For th is review , 
I will follow the evolutionary sequence outlined by Beuther et al.l 1 2007al l: 
Massive star-forming regions start as High-Mass Starless Cores (HMSCs), i.e., 
these are massive gas cores of the order a few 100 to a few 1000 Mq without 
any embedded protostars yet. In the next stage we have High-Mass cores with 
embedded low- to intermediate-mass protostars below 8M0, which have not 
started hydrogen burning yet. During that evolutionary phase, their lumi- 
nosity should still be dominated by accretion luminosity. Following that, the 
so called High-Mass Protostellar Objects (HMPOs) are still massive gas cores 
but now they contain embedded massive protostars >8M0 that have started 
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hydrogen burning which soon dominates the total luminosity of the sources. 
Hot Molecular Cores (HMCs) and hypercompact Hll regions (HCHlls) are 
part of that class. The last evolutionary stage then contains the final stars 
that have stopped accreting. While most ultracompact Hll regions (UCHlls) 
are Hkely part of the latter group, some of them may still be in the accretion 
phase and could then hence still harbor HMPOs. 

Definitions of a massive protostar: Another debate in massive star for- 
mation centers around the exact definition of a "massive protostar". If one 
followed the low-mass definition which basically means that a protostars is an 
object that derives most of its luminosity from accretion, then "massive pro- 
tostars" should not exist or only during a very short period of time because as 
soon as they are termed "massive" (>8Mq), their luminosity is quickly dom- 
inated by hydrogen burning. In this scenario, during the ongoing formation 
processes, one would then n eed to talk about "accretin g stars". This approach 
is for example outHned by IZinnecker &: Yorkd 1 20071 ) . A different definitio n 
for "massive protostars" is advocated, e.g., recently by lBeuther et al. I 2007al l: 
In this picture, a protostar is defined in the sense that each massive object 
that is still in its accretion phase is called a "massive protostar", independent 
of the dominating source of luminosity. This definition follows more closely 
the usual terminology of "proto", meaning objects that are not finished yet. 

Observational challenges: Whatever physical or chemical processes we are 
interested in massive star formation, one faces severe observational challenges 
because of the clustered mode of massive star formation and the on average 
large distances of a few kiloparsec. Therefore, high spatial resolution is a 
prerequisite for any such study. Furthermore, the early stages of massive star 
formation are characterized by on average cold gas and dust temperatures 
which are best observed at (sub)mm wavelength. Hence, most observations 
presented in the following are based on (sub)mm interferometer observations 
of young massive star-forming regions at different evolutionary stages. 

The main body of this article is divided into four sections deahng first with 
the initial conditions that are present prior to or at the onset of massive star 
formation. The next two sections deal with our current knowledge about the 
properties of potential massive accretion disks and the fragmentation behavior 
of massive gas clumps and cores. The following section will then outline the 
status and future possibilities of astrochemical investigations. Finally, I try 
to sketch the directions where current and future answers to the questions 
raised in the Abstract may lead to. 



2 The earliest stages of massive star formation 

What are the initial conditions of massive star formation? Until a few years 
ago, to address this question observationally in a statistical sense was close to 
impossible because we had no means to identify large samples of sources prior 
to or at the onset of massive star formation. The situation has changed sig- 
nificantly since the advent of the space-based near- and mid-infrared missions 
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that surveyed the Galactic plane starting with ISO and MSX, and now con- 
ducted with much better sensitivity and spatial resolution by Spitzer. These 
missions have revealed more than 10^ Infrared Dark Clouds (IRDCs), which 
are cold molecular cl ouds that are identified as shadows against the Galac- 
tic background (e.g.. lEgan et al.i. 1998; .Carey e t al. 2000; SimQll_et_aL, 2(306) . 
These clouds are characterized by on average cold temperatures (~15K), 
large masses (a few 100 to a few 1000 Mf)) and average densities of the order 
10"* -10^ cm" 1 fe.g. jRathborne et al. 2005: Sr idharan et al . 2005; Pillai et al.1 



2006al ). Although these clouds appear as dark shadows, they may be starless 
but they can also harbor embedded forming protostars. In fact, a statistical 
analysis of the percentage of starless IRDCs versus IRDCs with embedded pro- 
tostars will be an important step to understand the time-scales important for 
the earliest evolutionary stages. Currently, the statistical dat abase of in depth 



IRDC studies is still insufficient for such an estimate (e.g., Rathborne et al 



l2005Ll2006tlPillai et al.ll2006bHBeuther fc Steinackeijl20070Motte et al.ll2007t l 

however, it is interesting to note that until now no starless IRDC has been 
unambiguously identified in the literature, all detailed studies revealed em- 
bedded star formation processes. To first order, this triggers the interpre- 
tation/speculation that the high-mass starless core phase is likely to be ex- 
tremely shortlived. Future investigations of larger samples will answer this 
question more thoroughly. 



Figure 1: Sample IRDCs from lSridharan et al.l lj2005l ): MSX A-band (8^m) images 
(black is bright) with 1.2 mm emission contours: The first two numbers refer to the 
corresponding IRAS source and the third number labels the mm sub-sources. The 
five-pointed stars mark cores lacking good 1.2 mm measurements. 

In a recent spectral line study of a sample of 43 IRDCs (Fig.fTll. feeuther fc SridharanI 



1 20071 ) detected SiO(2-l) emission from 18 sources. Assuming that SiO is 
produced solely through sputtering from dust grains, and that this sample 
is representative for IRDCs in general, it indicates that at least 40% of the 
IRDCs have ongoing outfiow activity. Since the non-detection of SiO does 
not imply no outfiow activity, this number is a lower limit, and even a higher 
percentage of sources may harbor already ongoing star formation. The range 
of observed SiO line- widths down to zero intensity varied between 2.2 and 
65kms~^. While inclination effects and embedded objects of different mass 
could account for some of the differe nces, such effects are unlikely causing 



the whole velocity spread. Therefore, iBeuther &: SridharanI l|2007l l speculate 



whether the varying SiO line-widths are also indicators of their evolutionary 
stage with the smallest li ne-width close after the onse t of star formation ac- 
tivity. In the same study, Beuther fc SridharanI 1 2007 1 observed CH3OH and 



CH3CN. While CH3CN was detected only toward six sources, CH3OH was 
found in approximately 40% of the sample. The derived column densities 
are low of the order 10"^° with respect to H2. These values are consistent 
with chemical models of the earliest evolutionary stages of high-mass star for- 
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mation fe.g.. iNomura fc Millaiil2004l l. and the CH3OH abunda nces compare 
well t o recently reported values for low-mass starless cores (e-g.- iTafalla et alJ 



Zooming into selected regions in more detail, we studied one particu- 
larly interesting IRDC at high angular resolution with the Plateau de Bure 
Interferomet er and the Sp i tzer S p ace Telescope (IRDC 18223 -3, see Fig. [T] 
right panel; iBeuther et alJ l2005al : iBeuther fc Steinackerl 120071 ) . Combining 
the Spitzer mid-infrared data between 3 and 8/Ltm with the 3.2 mm long- 
wavelengths observations from the Plateau de Bure Interferometer (PdBI), 
we did not find any mi d-infrared c ounterpart to the massive gas core de- 
tected at 3.2 mm (Fig. [2| [Beuther et al. 2005al. However, we did detect three 
faint 4.5 features at the edge of the central 3.2 mm continuum core. Since 
emission features that occur only in the 4.5 /im but no other Spitzer band 
are usually attributed to sh ocked H2 emission from molecular outflows (e.g.. 



Noriega-Crespo et al.ll2004l ). we concluded that the region Hkely hosts a very 



young protostar that drives a molecular outflow but that is still too deeply 
embedded by to be detected by the Spitzer IRAC bands. This interpretation 
found further support by line-wing emission in older CO and CS data. Based 
on the inferred central source, we predicted that the region should have a 
strongly rising spectral energy distribution (SED) and hence be detected at 
longer wavelengths. As soon as the MIPSGAL mid- to far-infrared survey 
with Spitzer became availa ble, we the n could identify the central source at 24 
and 70 fim (Fig. [2^ .Beuther fc Steinack er 20071) . Combing the available mid- 
/far- infrared data with the long- wavelengths observations in the mm regime, 
it is possible to flt the SED with a two component model: one cold component 
(~15K and ~576Mq) that contains most of the mass and luminosity, and 
one warmer component (~51K and ^O.OIM©) to explain the 24 ^m data. 
The integrated luminosity of ~177L0 can be used to constrain additional pa- 
rameters of the embedded p rotostar from the tu rbulent core accretion model 
for massive star formation ( McKee fc Tan 2003l l. Following the simulations 
by iKrumholz et all ()2Q07bl ). the data of IRDC 18223-3 are consistent with a 
massive gas core harboring a low-mass protostellar seed of still less than half 
a solar mass with high accretion rates of the order lO~^M0yr~^ and an age 
below 1000 yrs. In the framework of this model, the embedded protostar is 
destined to become a massive star at the end of its formation processes. While 
this interpretation is attractive, it is not unambiguous, and especially the de- 
rived time-scale from this model appears short when comparing with recent 
outflow data that will be presented in the following section ( §3.ip . 



In summary, these observations indicate that the physical and chemical 
conditions at the onset of low- and high-mass star formation do not differ 
signiflcantly (except for largely different initial cloud clump masses and ac- 
cretion rates) , and that the time-scale for massive bound gas clumps to remain 
starless is likely relatively short. 
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Figure 2: IRDC 18223-3 images at different wavelengths from lSeuther fc Steinackerl 
The color scales show Spitzer images at various wavelength, and the 
contours show the 93 GHz (3.2mm) continuum emission observed with the PdBI 
ijSeuther et al.|[2005al l. The left panel presents a three-color composit e with blue 
3.6 fj,m, green 4.5 fim and red 8.0 fj,m (adapted from lSeuther et al.ll2005al l. The inlay 
zooms into the central core region. The middle and right panel show the Spitzer 
24 and 70 ^m images, respectively. The circles in each panel present the Spitzer 
beam sizes and the ellipse in the left panel presents the PdBI 3.2 mm continuum 
synthesized beam. A size-ruler is also shown in the left panel. 



3 Massive accretion disks? 

As mentioned in the Introduction, molecular outflows and accretion disks can 
be used to discriminate between the diflFerent formation scenarios for massive 
stars. Massive outflows have been subject to intense studies for more than a 



decade (e.g..iShepherd fc Chu rchwell 1996; Beuther et al. 2002b; Z hang et al.. 
20051 : Arce et all 20071 ) . and although there is considerable discussion about 



the details, we find a growing consensus that massive molecular outfiows are 
ubiquitous in high-mass star formation, and that collimated jet-like outflows 
do exist for m assive sources as well, at least during the very early evolu- 
tionary stages I Beuther fc ShepherdI 2005 ). The collimation of the outflows is 



Hkely to widen with ongoing evolution. Nevertheless, these data are consistent 
with the turbulent core model for massive star formation, whereas they are 
less easy to reconcile with the competitive accretion model because the latter 
is so dynamic that collimated structures likely could not survive very long. 
Furthermore, the existence of collimated outflows can only be explained by 
magneto-hydrodynamic acceleration of the jet from an underlying accretion 
disk. Hence, there is ample indirect evidence for massive accretion disks, how- 
ever, the physical characterization of d isks in massive star formation is still 



lacking largely the observational basis ijCesaroni et al.l 120071 ). The two main 
reasons for this are that the expected massive accretion disks are still deeply 
embedded within their natal cores complicating the differentiation of the disk 
emission from the ambient core, and that the clustered mode of massive star 
formation combined with the large average distances of the targets makes spa- 
tially resolving structures of the order 1000 AU a difficult observational task. 
In spite of these difficulties, the advent of broad spectral bandpasses allow- 
ing us to study several spectral Hnes simultaneously, as well as the improved 
spatial resolution of existing and forthcoming interferometers have increased 
the number of disk s tudies over the last few years. For a recent review see 



Cesaroni et all l|2007l ). 

Here, I will show three different examples of disk and/or rotation candi- 
dates in an evolutionary sense: It starts with a rotation and outffow investiga- 
tion of the previously discussed IRDC 18223-3 (|2|), then I will present recent 
data from the high-mass disk candidate in the HMPO IRAS 18089-1732, and 
finally observations from a massive disk candidate at a more evolved evolu- 
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tionary stage will be discussed. 



3.1 Rotation in IRDCs: the case of IRDC 18223-3 

As a follow-up of the Infrared Dark Cloud study of IRDC 18233-3 discussed in 
HI Fallscheer et al. (in prep.) observed the same region with the Submillimeter 
Array (SMA) in several spectral setups around 230 and 280 GHz covering 
outflow as well as dense gas tracers. Figure [3] shows a compilation of the 
CO (2-1) data, one CH3OH line and the dust continuum emission. The blue- 
and red-shifted C0(2-l) emission clearly identifies at least one large-scale 
outfiow in the north-west south-east direction. This is consistent with two 
of the 4.5 /zm emission features at the edge of the core (Fig. [21 left panel 
and inlay). There is another coUimated red-shifted feature to the south-west 
corresponding to the third 4.5 /zm feature, however, we do not identify a blue 
counterpart and refrain from further interpretation of that feature. Since we 
find for the main north-west south-east outfiow blue- and red-shifted emission 
on both sides of the continuum peak, the orient ation of the outfiow shou ld be 
close to the plane of the sky (see, e.g., models by lCabrit fc Bertout 1990f ). and 



hence the assumed underlying perpendicular rotating structure close to edge- 
on. Following the approach outlined in .Beuther et a l, (200_2b), the outfiow 
mass and outfiow rate are 13 and 3.5 x 10~^ M0yr~^, respectively. With 
the above derived core mass of ~576Mq (12]), this source fits well into the 
correlatio n between outfiow rat e and core mass previously derived for HMPOs 



(Fig. 7 in lBeuther et al.ll2002bf ) 



Fi gure 3i SMA observations toward IRDC 18223-3 (Fallscheer et al. in prep.). The 
left panel shows the blue- and red-shifted CO(2-l) emission as solid and dashed 
contours overlaid on the grey-scale 1.3 mm dust continuum emission. The central 
core is the same source as in Fig. [2] The right panel presents in grey-scale a 1st 
moment map (intensity weighted velocity distribution) of CH3OH overlaid with the 
1.1mm continuum emission. The empty and full circle are the synthesized beams 
of the line and continuum emission, respectively. 



While the outfiow rate is consistent with the accretion rate previously de- 
rived from the SED (|2|), discrepancies arise with respect to the a ge of the sys- 



tem. Although dynamical timescales are highly uncertain (e.g.. iParker et al 



199 It ) , the size of the molecular outfiow combined with a low inclination an- 



gle allows for at least a timescale estimate for the outfiow of the order a few 
10* yrs, well in excess of the value < lO'^ yrs previously derived from the SED 
appHed to models (|2]). Notwithstanding the large errors between the differ- 
ent estimates, the discrepancy of more than an order of magnitude appears 
real. How can we explain that? There is no clear answer to that yet, but 
a possibility is that the orientation of the disk-outfiow system with the disk 
close to edge-on absorbs a large amount of fiux distorting the SED on the 
Wien-side. If that were the case, the SED-estimated age could underestimate 
the age of the system. Another possibility to solve the discrepancy is that 
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the initial start of high-mass star formation may proceed slower, i.e., the first 
low-mass protostar(s) (destined to become massive or not?) form within the 
massive cores, and they already start driving outflows, but at that stage it 
is impossible to detect them in the near- to far-infrared because of the large 
extinction. In this picture at some point the high-mass star formation process 
would need to accelerate because otherwise t he massive stars can not form in 
the short time-scales of a few 10^ yrs (e.g., McKee fc TanI 2002f ). It is not 
clear why the whole process should start slow and what could trigger such 
acceleration later-on. Obviously, more theoretical and observational work is 
required to explain the different time-scales in more detail. 



Figure [3] (right panel) zooms into the central core and shows dust contin- 
uum emission as well as the velocity structure of the dense central gas observed 
in CH30H(6o,i — 5o,i) with a lower level excitation level of E'low — 34.8 K. In- 
terestingly, both the continuum and the spectral line emission are elongated in 
the north-east south-west direction perpendicular to the main molecular out- 
flow. While the continuum emission shows three resolved emission features, 
CH3OH exhibits a smooth velocity gradient across the source spanning ap- 
proximately 3kms~^. The CH3OH line-width FWHM toward the continuum 
peak is 2.1kms~^. The blue-redshifted features in the north-west are likely 
part of the molecular outflow and one sees even a slight elongation of the con- 
tinuum emission in that direction. While CH3OH is a well-known shock tracer 



and hence regularly found within molecular outflows (e.g., iBachiller et al 



20011 1 , it is more of a surprise that we flnd it in an elongated structure likely 



associated with rotation and infall perpendicular to the outflow. The extent 
of this structure is large with ~6.5" corresponding to more than 20000 AU at 
the given distance of 3.7 kpc. Although we have no methanol isotopologue in 
the setup to exactly determine the opacity of the line, a low-energy transition 
like this is likely to be optically thick. Hence, we are tracing some of the outer 
rotating structures, probably correspondin g to a larger sca le rotating and po- 
tentially infalling/inspiralling toroid (e.g.. ICesaroni et al.li2005,: KetQ,2007i) . 
The small velocity spread across the structure as well as the relatively narrow 
CH3OH Hue- width toward the core center are also consistent with tracing 
outer structures because due to momentum conservation rotating structures 
should have lower velocities further out. 



Notwithstanding that we do not exactly know the age of IRDC 18223-3, 
its non-detection up to 8^m puts it at an early evolutionary phase prior to 
the better studied HMPOs and Hot Molecular Cores. Our data clearly show 
that even at such early stages molecular outflows and rotating structures 
perpendicular to that have been developed, and it is likely that closer toward 
the core center, one will flnd a real accretion disk. To investigate the latter in 
more detail, higher angular resolution observations of an optically thin dense 
gas tracer are required. 



8 



3.2 The HMPO disk candidate IRAS 18089-1732 



As a more evolved massive disk candidate, we have studied intensely over 
the last few years the HMPO IRAS 18089-1732. This source at a distance 
of 3.6 kpc with a luminosity of 10^'^ L© is part of a large sample of HM- 
POs, it hosts H2O and Class II CH3OH maser and has s trong molecular 



line emission indicative of e n embedded Hot Molecular Core (ISridharan et al 
2QO2I: iBeuther et al.ll2002al l. During early SMA observations ISeuther et al 



ll2n04d lal \20(i^ identified in SiO a molecular outfiow in the north-south di- 



rection, and perpendicular to that in HCOOCH3 a velocity gradient on scales 
of a few 1000 AU. Although these data were indicative of rotation and an 
underlying massive accretion disk, the observations did not allow us to char- 
acterize the structure in more detail because of a lack of spatial resolution. 
Therefore, we observed IRAS 18089-1732 now in high-energy transitions of 
NH3 at 1.2 cr n wavelength with the V LA and the ATCA at a spatial resolu- 
tion of 0.4" llBeuther fc Walshll2008l l. These NH3(4,4) and (5,5) lines have 
a two-fold advantage: Their high excitation levels (> 200 K) ensure that we 
are tracing the warm inner regions and are less confused by the surrounding 
cold envelope, whereas the cm wavelengths regime is less affected by high 
optical depth of the dust emission in high column density regions and may 
hence be particularly well suited for massive disk studies (e.g. Krumholz et al.l 



2007al) . Figure [4] presents an integrated image and a 1st moment map (inten- 



sity weighted velocity) of the corresponding VLA observations. 



Figure 4: T he left pane l shows the VLA NH3(5,5) emission integrated from 31 
to 37kms"^ ijBeuther fc W alsh 20oi l. The right panel presents the corresponding 
1st moment map contoured from 31.5 to 36.5 kms~^ (step lkms~^). The white- 
black dashed contours show the 1.2 cm c ontinuum emission. T he asterisks mark the 
position of the submm continuum peak ijBeuther et al.|[2005d l. and the synthesized 
beams are shown at the bottom- left (grey NH3 and dashed 1.2 cm emission). 

The 1st moment map confirms the previously assessed velocity gradient 
in east-west direction perpendicular to the molecular outflow. The NH3 Hne- 
width FWHM toward the central core is 4.7kms~^, signiflcantly broader than 
that of IRDC 18223-3 f ^3.ip . In the simple picture of equilibrium between 
gravitational and centrifugal force, the rotationally bound mass would be 
~37M0, of the same order as the whole gas mass as well as the mass of 
the central source (of the order 15 M©). Furthermore, the position-velocity 
diagram is not consistent with Keplerian rotation. It even shows indications 
of super-Keplerian motion, which is expected for very massive disks where 
the rotation profile is not only de termined by the mass o f the central object 



but also by the disk itself (e.g., iKrumholz et al.M2007bf ). Hence, the new 



VLA and ATCA data clearly confirm the previous assessment of rotation 
perpendicular to the outfiow /jet, however, the kinematic signatures of that 
rotating structure are not consistent with a Keplerian disk Hke in low-mass 
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star formation, but they show additional features which can be produced 
by massive self-gravitating disks as well as by infalling gas that may settle 
eventually on the disk. In addition to this, the detection of the high-excitation 
lines in the rotating material indicates high average gas temperatures >100 K 
for the disk-like structures, well in excess of typical g as temperatures in low- 



mass disks of the order 30 K (e.g., iPietu et al.l 120071 ). Moreover, we detect 



double-lobe cm continuum emission close to the core center where the two 
lobes are oriented in north-south direction parallel to the outflow identifled in 
SiO. With respect to previous data at longer wavelength, we flnd a spectral 
index at cm wavelength of 1.9, consistent with an optically thick jet (|Reynolds 

It will be interesting to further zoom into the innermost regions with 
future instruments like ALMA and eVLA to asses whether the quantitative 
deviations from typical low-mass accretion disks continue down to the smallest 
scales, or whether we will find Keplerian disk structures as known from their 
low-mass counterparts. 



3.3 A more evolved massive disk candidate? 

Moving along in the evolutionary sequence, we have recently identified a po- 
tential disk around a more evolved candidate young stellar object (Quanz et 
al. in prep.). The source labeled so far mdcl (massive disk candidate) was 
identified serendipitously during a near-infrared wide-field imaging campaign 
on Calar Alto via its K-band cone-like nebulosity and a central dark lane 
(Fig. HI). First single-dish bolometer and spectral line measurements revealed 
a 1.2 mm flux of 12 mJy and a velocity of rest of '^51 km s~^. The latter value 
indicates a kinematic distance of ^--^5 kpc, consistent with distances of a few 
UCHii regions in the surrounding neighborhood. To investigate this object 
in more detail, we recently observed it with the SMA at 1.3 mm wavelength 
mainly in the mm continuum and the ^^CO/C^^O spectral line emission. Fig- 
ure [5] presents an overlay of the SMA data with the K-band nebulosity, and 
a few points need to be stressed: 



Figure 5: The grey-scale in both panels shows the K-band near-infrared nebulosity 
observed for this massive evolved disk candidate (Quanz et al. in prep.). The con- 
tours are the corresponding SMA mm observations where the left panel shows the 
1.3 mm continuum emission, and the right panel in black and white contours the 
blue-shifted ^^C0(2-l) and the integrated C^**0(2-l) emission, respectively. 

(a) Although spatially unresolved with a synthesized beam of ^ 4.0" the 
1.3 mm continuum peak exactly coincides with the infrared dark lane consis- 
tent with the large column densities of the proposed disk-like structure. The 
flux measured with the SMA is 12 mJy like the previous single-dish measure- 
ments. This indicates that we have no surrounding dust /gas envelope but 
rather an isolated central structure. Assuming optically thin dust emission 
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at 50 K, the approximate gas mass of the central structure is ^5Mq. (b) We 
detect blue-shifted CO (2-1) spatially well correlated with the K-band nebu- 
losity north of the dark lane. This confirms the initial interpretation of that 
feature to be due to an outfiow. (c) The integrated C^^0(2-l) emission is 
elongated perpendicular to the outflow observed in CO and K-band contin- 
uum emission. The line-width FWHM of the C^*0 emission is narrow with 
'^0.8kms~^, however, the spatial extent of this structure is large, of the order 
2 X 10-* AU. 

While the low gas mass and the missing more massive gas envelope could be 
interpreted in the framework of a low-mass source, such large disk-structures 
as indicated by the C^^O e mission are not known from typical low-mass disk 



sources (e.g.. ISimon et al.ll200ff ). Therefore, these observations can also be 



interpreted as a remnant disk/torus around an intermediate to high-mass 
(proto)star that has already dispersed much of its envelope. Although C^^O 
is detected only in two channels, these show a clear velocity shift, and the 
small line-width may be due to the lower rotational motions on large scales 
assuming momentum conservation in rotating, potentially Keplerian struc- 
tures. 



Synthesizing the three example sources shown here, it is interesting to 
note that the line-widths are small in the youngest and the supposed to be 
oldest source, whereas they are large in the HMPO which should be in its 
main accretion phase. In an evolutionary picture this can be interpreted that 
at early evolutionary stages infall, turbulence and rotation are not yet that 
vigorous. Then in the main accretion phase, infall, rotation and outflow pro- 
cesses strongly increase the line-width. And flnally, when the accretion stops, 
the envelope and disk slowly disperse and one observes only a remnant struc- 
ture with small line-width in the outer regions. This scenario is speculative, 
however, the number of disk candidates is steadily increasing, and since we 
start sampling more evolutionary stages, we are getting the chance to address 
disk evolution questions in high-mass star formation as well. 



4 Fragmentation in high-mass star formation 



How massive gas clumps fragment is one of the key questions if one wants 
to understand the formation of the Initial Mass Function, and as outlined 
in m the two main massive star formation theories predict differences in the 
early fragmentation processes. In the following I will present several examples 
of fragmenting massive cores addressing issues about fragmentation on the 
cluster-scale, fragmentation of smaller groups, potential proto-trapezia, and 
the determination of density structures of sub-sources within evolving clusters. 
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4.1 Resolving the massive proto-cluster IRAS 19410+2336 



To address fragmentation processes at early evolutionary stages high angu- 
lar resolution at (sub) mm wavelengths is t he to ol of choice to resolve the 
relevant substructures. iBeutherfc Schilkd tOOi \ resolved the young mas- 
sive star-forming region IRAS 19410+2336 (distance ~2.1 kpc and luminosity 
^IO'^Lq) at 1.3 mm wavelength with the PdBI at approximately 2000 AU lin- 
ear resolution into 24 sub-sources. Although from a statistical point of view 
such numbers cannot compete with the clusters exceeding 100 or even 1000 
stars observed at optical and near-infrared wavelength, this is still one of the 
prime examples of a spatially resolved massive proto-cluster. Assuming that 
all emission features are due to cold dust emission from embedded protostars, 
they were able to derive a core-mass function. With a power-law slope of 
-2.5, this core mass function is consistent with the Salpet er IMF slope of - 
2.35 ()Salpeted[l955h . Therefore. [Beuther fc Schilkd ()2004h interpreted these 
observations as support fo r the turbulent fragmentation put forth by, e.g., 
Padoan fc Nordlundl l|2002f l. 



A few caveats need to be kept in mind: While iBeuther &: Schilkd ll2004h 
assumed a uniform gas temperature for all sub-sources, it is more likely that 
the central peaks are warmer than those further outside. This issue can be 
addressed by spectral line emission with temperature sensitive molecules (e.g., 
H2CO) which is an ongoing project by Rodon et al. (in prep.). Furthermore, 
the assumption that all mm continuu m peaks are o f pro- or pre-stellar na- 
ture is not necessarily always valid, e.g., iGueth et al] ()2003l ) orlBeuther et al 



I 2004bl l have shown that mm continuum emission can partly also be caused 
by coUimated jets. However, only the central source is detected at cm wave- 
length and collimated jets should be detectable at cm wavelengths as well. 
Therefore, we believe that jets should not affect the analysis much. 



Independent of the caveats, it is surprising that IRAS 19410+2336 is still 
the only young massive star-forming region that is resolved in >10 sub-sources 
in the mm continuum emission. While this can be explained to some de- 
gree by the exceptionally good uv-coverage obtained for the given observa- 
tions, which results in a good sampling of spatial structures, we also need to 
consider whether different modes of fragmentation may exist. Similar high- 
spatial-resolution studies of more proto-clusters spanning a broad range of 
luminosities are required to tackle this question in more detail. Another in- 
teresting question is associated with the spatial filtering of interferometers and 
the corresponding large-scale emission: Many interferometric (sub)mm con- 
tinuum studies of massive star-forming regions filter out of the order 90% of 
the flux, hence, large amounts of the gas are distributed on larger scales, usu- 
ally > 10". The question remains whether this gas will eventually participate 
in the star formation process or not? 
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4.2 Fragmentation of potential proto-trapezia 



4.2.1 The enigmatic proto-trapezium W3-IRS5 



The W3-IRS5 region is one of the prototypical high-mass star-forming re- 
gions with ^ 10^ Lq at a distance of ~1.8kpc that shows fragmentation on 
scales of the order 100 AU observed at near-infrared as well as cm wave- 

2nni 



lengths (jMegeath et al 



van der Tak et al.ll2005f ). However, not much 



was known about the cold dust and gas emission. Therefore, we observed the 
region with the PdBI at 1.3 and 3.5 mm wavelengths with the new extended 
basehnes resulting in an unprecedented spatial resolution of ~ 0.37" (Rodon 
et al. in prep.). Figure [6] shows a compilation of the 1.3mm continuum data 
and the SiO(5-4) and 802(222, 20 22i 21) spectral line emission. 



Figure 6: PdBI observations of tlie W3-IRS5 system from Rodon et al. (in prep.). 
Tlie left panel sliow s the 1.3 mm c o ntinu um emission, and tlie stars mark near- 
infrared sources from lMegeath et al.l l|2005l ). The middle panel presents as solid and 
dotted contours the blue- and red-shifted SiO(5-4) emission overlaid on the grey- 
scale 1.3 mm continuum emission. The right panel finally shows the 1st moment 
map of 802(222,20 22i,2i in grey-scale with the 1.3 mm continuum contours. 

The mm continuum emission resolves the W3-IRS5 region into five sub- 
sources, where four of them are coincident with near-infrared and cm emission 
peaks. Three of the sources are clustered in a very small projected volume 
of only ^^2000 AU. With this high spatial resolution we find extremely large 
average column densities of the order a few times 10^"* cm~^ which corre- 
sponds to visual extinctions Ay between 5 x lO'^ and 10^ averaged over the 
beam size. Such extinctions should be far too large to allow any detection 
at near - infrared wavelengths , nevertheless, near-infrared counterparts are de- 



tected l|Megeath et al.M2005f ). This conundrum can likely be explained by 
the detection of several SiO outfiows in the field. In particular, we find very 
compact blue- and red-shifted SiO emission toward the two main mm peaks, 
where the blue- and red-shifted emission is barely spatially separated (Fig. [6] 
middle panel) . Since the overall time-scale of th e W3-IRS5 outflow sy stem is 
relatively large (of the order a few times 10^ yrs.lRi dge fc Moorell200ll ). these 
compact features are unlikely from very young outflows, but they indicate 
that the outflows are oriented almost along the line of sight. The opening 
cones of the outflows are the likely cause that emission from close to the pro- 
tostars can escape the region and hence make them detectable at near-infrared 
wavelengths. 

The right panel of Fig. [6] shows the 1st moment map of SO2 (intensity 
weighted velocities) which encompasses the mm continuum peaks. The co- 
herent velocity fleld over the sub-sources is a strong indicator that the system 
is a boun d structure and not some unbound chance alignment within the 



field (e.g., Launhardt 2Qq3). In addition to the general velocity gradient from 



the south-east to the north-west, one tentatively identifies velocity gradients 
across the two strongest mm continuum peaks. Since we do not know the 
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exact orientation of the outflows with respect to the SO2 rotation axis, it 
is not yet possible to identify these structures with disk-Hke components as 
in 131 Future observations in different tracers may help to shed more light 
on the rotational structure associated with each sub-source. It should also 
be noted that the line-width FWHM toward the mm continuum peaks varies 
between 6.2 and 7kms~^ larger than those reported in |3l While the larger 
line-width compared with the IRDC and the more evolved source may be 
explained by the evolutionary sequence sketched at the end of ^ the larger 
FWHM compared to the HMPO may have different reasons, among them 
are the larger luminosity of W3-IRS5, its multipHcity compared with the so 
far unresolved source IRAS 18089-1732, and also the molecular species, be- 
cause SO2 should be more affected by shocks than the NH3 line used for the 
IRAS 18089-1732 study. In addition to this, the SO2 moment map exhibits a 
velocity discontinuity with a velocity jump of the order 4kms~^ south-east of 
the mm continuum peaks. What is the cause of this discontinuity, is it asso- 
ciated with the original core formation and a shock within converging flows, 
or is it of different origin? 

In summary, the combination of high-spatial-resolution observations of the 
continuum emission in addition to outflow and dense gas tracers allows us to 
characterize many physical properties of this proto-trapezium system with 
respect to its multiple components and their outflow and rotation properties. 



4.2.2 Fragmentation of the hot core G29 



The hot core G29.96 located right next to a well-known cometary Hii re- 
gion comprises another example of several protostellar submm continuum 
sources within the innermost center of a high-mass star forming region (dis- 
tance '^6kpc, luminosity 9 x 10^ L©). High-spatial resolution observations 
with the SMA in its most extended conflguration yielded a spatial resolu- 
tion of 0.36" X 0.25" in the submm conti nuum at ^^348 GHz, corresponding 



to a linear resolution of 2000 AU (Fig. [3 iBeuther et al.ll2007d . the Hne data 
will be discussed in ^5.ip . The Hot Molecular Core previously identified in a 
high-excitation NH3 line ( Cesaroni et al.|[l998h is resolved by these new data 
into four sub-sources within a projected diamet er of ^6900 AU. Assu ming 



that the emission peaks are of protostellar nature. iBeuther et al.1 (|2007ci ) esti- 
mated a protostellar density of '--^ 2 x 10^ protostars/pc~^. This is considered 
a lower limit since we are limited by spatial resolution, sensitivity and pro- 
jection effects. Nevertheless, such a protostellar density is about an order of 
magnitude higher t han values usually reported for star-forming regions (e.g., 
Lada fc Ladal I2OO3I ) . Although this value is still about an order of magni- 



tude lower than protostellar de nsities that would be required in the niergin g 
scenario for massive stars (e.g., Bonnell et al. 2004; Ballv fc Zinnecker 2005f l. 
it is interesting to note that increasingly higher protostellar densities are re- 
ported when going t o younger sources and better angular resolution (see also 



Megeath et al. 20051 ). This allows us to speculate whether future observa- 



tions with better spatial resolution and sensitivity toward extremely massive 
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star-forming regions will reveal protostellar densities that may be sufficient 
to make mergers possible. While such a detection would not be a proof for 
mergers to exist, it will certainly be important to verify whether the required 
initial conditions do exist at all. 



Figure 7: Compilat ion of data toward the UCHii/hot core region G29.96 from 
iBeuther et"aD l2007d . The dashed contours present the cometary UCHll regions 
whereas the full con tours show the older NH3 observation from the hot core 
ijCesaroni et al.l[l994l l. The grey-scale with contours then present the new high- 
resolution (0.36" X 0.25") submm continuum data from the SMA. 



4.3 Density structure of sub-sources - IRAS 05358+3543 



As a final example for the potential of (sub)mm continuum studies, I present 
the recent multi-wavelength investigation of the HMPO IRAS 05358+3543. 
This region at a distance of 1.8 kpc with a luminosity of 10'^-^ Lq was observed 
in a combined effort with the PdBI and the SMA at arcsecond resolution in 



four wavelength bands (3. 1 and 1.2 mm, and 875 and 438 ^m. IBeuther et al 



iLeurim et al.ll2007l ). While many details about the sub-structure of 
the forming cluster can be derived, her e, I will discuss only tw o results. 

Based on the multi- wavelength data. IBeuther et ~al\ ( 2007bh fitted the spec- 
tral energy distribution on the Rayleigh- .Jeans side of the spectrum (Fig. [HI). 
While the main source can well be fitted by a typical protostellar spectrum 
consisting of free-free emission at long wavelength and a steep fiux increase 
at shorter wavelength due to the dust emission, another sub-source did not 
fit at all into that picture. In particular the shortest wavelength data-point 
at 438 /im shows significantly lower fiuxes than expected for a typical proto- 
star. The most likely explanation for this effect is that we are deahng with 
a very cold source and that therefore we are already approaching the peak of 
the spectral energy distribution. The data allowed us to estimate an upper 
limit for the dust temperature of < 20 K. Since we are also not detecting any 
other line ernission from this core (mainly from typical hot core molecules. 



iLeurini et al.ll2007l l. it may well be a starless core right in the vicinity of an 
already more evolved massive protostar. Further investigations of this sub- 
source in typical cold gas tracers Hke N2H+ or NH3 are required to test this 
proposal. Independent of whether this source harbors an embedded proto- 
star or not, these observations show the importance of short wavelength data 
at high spatial resolution if one wants to differentiate between critical core 
parameters like the dust temperature. 

Another physical parameter which has so far not been observationally con- 
strained for massive star formation, is the density profile of individual sub- 
sources. While de nsity profiles of low-mass star-forming cores have well been 
characterized (e.g., Motte et al.lll998t Ward-Thompson et al.lll999l : Andre et al 
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Figure 8: The le ft panel presents th e SED toward the coldest sub-source in 
IRAS 05358+3543 ijSeuther et al.ll2007bl l. The parameters of the fits are marked 
in the figure. The right panel shows intensities averaged in uv-annuli and plotted 
versus the baseline-length for different sub-sources and wavelengths. Most can be 
well fitted by power-law distributions. 



2000I 1. in high-mass star formation, density profiles were until now only de- 



rived with single-dish observatio n covering scales of t h e whole cluster but 



not individual sub-sources (e .g., iBeuther et al.l 2002at Mueller et al.l l2002t 



Hatchell fc van der This is partly due to the technical problem 



of interferometer observations that filter out large amounts of the flux and 
hence make density profile determi nations from their ima ges extremely unre- 
liable. To overcome this problem, iBeuther et al.l (|2007a ) analyzed the data 
directly in the uv-domain prior to any fourier transformation. Figure [8] shows 
the corresponding plots of the observed intensities versus the uv-distance for 
three sub-sources in three wavelengths bands, respectively. The observations 
cannot be fitted with Gaussian distributions, but much better fits are achieved 
with power-law distributions. These power-laws in the uv-domain can directly 
be converted to the corresponding power-laws of the intensity profiles in the 
image plane. Assuming furthermore a temperature distribution T oc r~'^-'^ 
we can now infer the density profiles of individual sub-sources of the evolv- 
ing cluster. The derived density profiles p oc r^P have power-law indices p 
between 1.5 and 2. Although this result is similar to the density profiles pre- 
viously determined for low-mass cores, to our knowledge this is the first time 
that they have been observationally constrained for resolved sub-sources in 
a massive star-forming region. The density str ucture is an importan t input 
parameter for any model of star formation (e.g., McKee &: Tan 2003l l. 



5 Astrochemistry 

5.1 Toward a chemical evolutionary sequence 

Astrochemistry is a continuously growing field in astronomy. Although Hne- 
survey style studi e s of different source s have existed for quite some time (e.g., 
Blake et al.lll987l : ISchilke et al these studies had usually been per- 



formed with single-dish instruments averaging the chemical properties over 
the whole cluster-forming regions. Since the advent of broadband receivers 
at interferometers like the SMA, it is now also possible to perform imaging 
spectral line surveys that allow us to spatially differentiate which molecules 
are present in which part of the targeted regions, for example, the spatial 
differ entiation between n itrogen- and oxygen-b earing molecules in Orion-KL 
(e.g., iBlake et al. 19961 : Beuther et al. 2005b[ ). In addition to the spatial 



analysis of individual regions, we are also interested in analyzing how the 
chemistry evolves in an evolutionary sense. As an early step in this di- 
rection we synthesized SMA observations that were observed in the same 
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spectral setup around 862 toward four massive star-forming regions over 
the last few years (Beuther et al. subm.). These four regions comprise a 
range of luminosities between lO^'^L© and 10^ L©, and they cover different 
evolutionary stages from young High-Mass Protostellar Objects (HMPOs) 
to typical Hot Molecular Cores (H MCs): Orion-KL: HMC, L ~ IO^Lq, 
D - 45kDC llBeuther et al.l [2005bl): G29.96: HMC, L - 9 x IO^Lq, D - 
6kpc iBeuther et al."2007d); IRAS 23151, HMPO, L ~ lO^L©, D - 5.7kpc 

HMPO, L ~ lO^-^L©, D - 1.8 kpc 
Smoothing all datasets to the 



Beuther et al. 


2007d^ 


Beuther et al. 


20n7bl: 



IRAS 05358: 



Leurini et all 120071 ) 



same linear spatial resolution of 5700 AU, we are now capable to start com- 
paring these different regions. Figure [9] presents typical spectra extracted 
toward the HMC G29.96 and the HMPO IRAS 23151. 




Rest Frequency (MHz) 



Figure 9: SMA spectra extracted toward two massive star-forming regions (G29.96 
top row & IRAS 23151+5912 bottom row, Beuther et al., subm.). The spectral 
resolution in all spectra is 2km/s. The left and right column show the lower and 
upper sideband data, respectively. 



A detailed comparison between the four sources is given in a forthcoming 
paper (Beuther et al. subm.), here we just outline a few differences in a qual- 
itative manner. 

• The HMCs show far more molecular lines than the HMPOs. Orion-KL and 
G29.96 appear similar indicating that the nature of the two sources is Hkely 
to be comparable as well. Regarding the two HMPOs, the higher luminos- 
ity one (IRAS 23151) shows still more Hnes than the lower-luminosity source 
(IRAS 05358). Since IRAS 05358 is approximately three times closer to us 
than IRAS 23151, this is not a sensitivity issue but it is likely due to the dif- 
ferent luminosity objects forming at the core centers. 

• The ground-state CH3OH lines are detected toward all four sources. How- 
ever, the vibrational-torsional excited CH3OH are only strongly detected to- 
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ward the HMCs Orion-KL and G29.96. Independent of the luminosity, the 
HMPOs exhibit only one CH3OH f t = 1 line, which can easily be explained 
by the lower average temperatures of the HMPOs. 

• A more subtle difference can be discerned by comparing the SO2 and 
the HN^^c/cjjgCjj^CN line blend near 348.35 GHz (in the upper sideband). 
While the SO2 line is found toward all four sources, the HN13C/CH3CH2CN 
line blend is strongly detected toward the HMCs, but it is not found toward 
the HMPOs. In the framework of warming up HMCs, this indicates that 
nitrogen-bearing molecules are either released from the grains only at higher 
temperatures, or they are daughter molecules which need some time during 
the warm-up phase to be produced in gas-phase chemical networks. In both 
cases, such molecules are not expected to be found much prior to the forma- 
tion of a detectable HMC. 

• Comparing the spatial distribution of different molecules, we find, e.g., that 
C'^''S is observed mainly at the core edges and not toward the submm con- 
tinuum peak positions. This difference can be explained by temperature- 
selective desorption and successive gas-phase chemistry reactions: CS desorbs 
early from the grains at temperatures of a few 10 K and should peak during 
the earliest evolutionary phases toward the main continuum sources. Sub- 
sequently when the core warms up to '^100 K, H2O desorbs and dissociates 
to OH. The OH then quickly reacts with the sulphur to form SO and SO2 
which should then peak toward the main continuum sources. This is what 
we observe in our data. The fact that the C'^'^S peaks are offset from the 
submm continuum condensations even toward the younger sources is due to 
their evolutionary stage where they have already heated up their central re- 
gions to more than 100 K. Even younger sources are required to confirm this 
scenario. 



5.2 C2H as a tracer of the earliest evolutionary stages? 

In an effort to study a larger source sample with respect to its chemical evo- 
lution, we observed 21 massive star-forming regions covering all evolutionary 
stages from IRDCs via HMPOs/hot cores to UCHll with the APEX tele- 
scope at submm wavelengths (Beuther et al. subm.). While most spectral 
fines were detected mainly toward the HMPO/hot core sources, the ethynyl 
molecule C2H is omni-present toward all regions. To get an idea about the 
spatial structure of ethynyl, we went back to an older SMA data-set targeting 
the HMP O IRAS 18089-1732 a t the same frequency around 349.4 GHz of the 



C2H line ijBeuther et al.ll2005d l. Because we were not able to image the spa- 
tial distribution of C2H at that time, we now restricted the data to only the 
compact configuration allowing us to better image the larger-scale distribution 
of the gas. Figure [10] presents the resulting molecular line map, and we find 
that C2H is distributed in a shell- like fashion around the central protostellar 
condensation. Comparing this with all other imaged molecules in the original 
paper, only C2H exhibits this behavior. To better understand this effect, we 
ran a set of chemical models in ID for a cloud of 1200 Mq, a density power-law 
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p (xr~^-^ and different temperature distributions T cx r"^. A snapshot of these 
models after an evolutionary time of 5 x 10** yrs is presented in Figure [TOl The 
models reproduce well the central C2H gap in IRAS 18089-1732 which should 
have approximately the same age. 



Figure 10: The left panel presents in grey-scale the C2H emission and in thick solid 
contours the corresponding submm continuum from the SMA toward the HMPO 
IRAS 18089-1732 (Beuther et al., subm.). The right panel shows a chemical model 
explaining the decreased emission toward the core center after approximately 5 x 
10* yrs. The parameter q denotes the temperature power-law index, and the T 
values refer to the temperature at the core edge or to isothermal values (g = 0). 



While these models reproduce the observations, they give also predictions 
how the C2H emission should look like at different evolutionary times. In par- 
ticular, C2H forms quickly early on, also at the core center. Since not many 
molecules exist which do not freeze out and are available to investigate the 
cold early phases of massive star formation (valuable exceptions are, e.g., NH3 
or N2H"'"), the detection of C2H toward the whole sample in combination with 
the chemical models triggers the prediction that C2H may well be an excellent 
molecule to investigate the physical conditions of (massive) star-forming re- 
gions at very early evolutionary stages. High-spatial-resolution observations 
of IRDCs are necessary to investigate this potentially powerful astrophysical 
tool in more detail. 



5.3 Employing molecules as astrophysical tools 

While the chemical evolution of massive star-forming regions is interesting in 
itself, one also wants to use the different characteristics of molecular lines to 
trace various physical processes. In contrast to molecules like SiO and CO 
that are well-known outflow /jet tracers, the task gets more difficult searching 
for suitable accretion disk tracers. Investigating our sample and disk claims 
in the literature, one finds that in many cases exclusively one or the other 
molecule allows the investigation of rotational motion, whereas most other 
molecular lines remain without clear signatures. For example, the HN^'^C 
line discussed above ( ^5.ip traces rotation in the hot core G29.96 but it is 
not even detectable in younger sources. The o ther way arouiid, C^ '^S traced 
disk rotation in the young HMPO IRAS 20126 (ICes aroni et~ani2005l l. but not 



anymore toward more evolved sources ( §5.ip . These differences imply that 
one will unlikely find a uniquely well suited molecular line allowing the study 
of large samples of massive accretion disks, but that one has to select for each 
source or source class the suitable molecule for detailed investigations. 

In the following, I give a short table with molecules and their potential 
usefulness to study different physical processes. This table (1) is restricted to 
molecules with spectral lines at cm/ (sub) mm wavelengths and does not claim 
any kind of completeness, it should just serve as a rough overview and it only 
Hsts the main isotopologues of each species. 
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OH Zeeman effect, magnetic fields, maser signpost of star 
formation 

CO General cloud structure, outflows 

SiO Shocks due to jets/outflows 

CO^ Far-UV radiation from embedded protostars 

CS Dense gas, rotation, also outflows 

CN Photodominated regions, Zeeman effect, magnetic flelds 

SO Shocks, dense gas 

H2O Shocks and hot cores, rotation (H2^0), maser signpost of 

star formation 

HDO Deuterium chemistry 

H2D+ Cold gas, pre-stellar cores, freeze out 

HCN Dense cores, also outflows 

HNC Dense cores, rotation (HN^^C) 

HCO"*" Outflows, infall, cosmic rays, ionization degree, 

dense gas (Hi3C0+) 

SO2 Shocks, dense gas 

C2H Early evolutionary stages ( §5.2p 

N2H+ Early evolutionary stages 

N2D"'" Deuteration, freeze out 

HaO"*" Cosmic rays 

H2CO Dense gas, temperatures 

NH3 Cold and hot cores, rotation, temperatures 

CH3OH Shocks, young rotating structures? ( ^3.ip . temperatures, 

maser signpost of massive star formation 

CH3CN Hot cores, temperatures, rotation 

CH3CCH Dense gas, temperatures 



HCOOCH3 Hot cores, rotation 
CH3CH3CN Hot Cores 



Table 1: A few useful molecules and some of their potential applications. 



6 Conclusions and summary 

This article tries to outline how far we can currently constrain physical and 
chemical properties in massive star formation using (sub)mm interferometry. 
Coming back to the original questions raised in the abstract: (a) What are 
the physical conditions at the onset of massive star formation? (b) What 
are the characteristics of potential massive accretion disks and what do they 
tell us about massive star formation in general? (c) How do massive clumps 
fragment, and what does it imply to high-mass star formation? (d) What 
do we learn from imaging spectral line surveys with respect to the chemistry 
itself as well as for utilizing molecules as tools for astrophysical investigations? 

Can we reasonably answer any of these questions with confldence? There 
are no clear-cut answers possible yet, however, the observations are paving a 
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way to shedding light on many of the issues, and one can try to give tentative 
early answers. The following is a rough attempt to outline the directions for 
current and future answers in these fields: 

(a) Massive gas clumps prior or at the onset of high-mass star formation 
are characterized by cold temperatures of the order 15 K and small line- widths 
indicative of a low level of turbulence. Their molecular abundances appear 
comparable to those of low-mass starless cores. Interestingly, the outflow 
detection rates toward IRDCs are high, and no genuine High-Mass Starless 
Cores have been reported in the literature yet. Although the statistical basis 
is not solid enough yet, this ahows us to speculate that the high-mass starless 
phase is likely to be very shortlived. 

(b) The detection of a real accretion disk around a massive protostar still 
remains an open issue. However, we find many rotating structures in the 
vicinity of young massive star-forming regions all the way from IRDCs to Hot 
Molecular Cores. These structures are on average large with sizes between 
1x10-^ and 2x10'' AU, and they ha\ 'c masses of the order of the central 
protostar. Hence, most of them are not Keplerian accretion disks but rather 
some larger-scale rotating/infalhng structures or toroids that may feed more 
genuine accretion disks in the so far unresolved centers of these regions. 

(c) Fragmentation of massive star-forming regions is frequently observed, 
and the core mass function of one young region is consistent with the Initial 
Mass Function. However, caveats of unknown temperature distributions or 
missing flux on larger scales may still afl^ect the results. Furthermore, we 
flnd proto-trapezium like structures which show multiple bound sources on 
small scales of a few 1000 AU implying protostellar densities of the order 10^ 
protostars/pc^'^. Such densities are stih not sufficient to allow coalescence, 
however, it may be possible to find even higher protostellar densities with the 
improved observational capabilities of future instruments. Although mergers 
do not appear necessarjf to form massive stars in general, they still remain a 
possibility for the most massive objects. 

(d) Astro-chemistry is a young branch in astrophysical research, and we 
arc currently only touching the surface of its potential. The different paths to 
follow in the coming years are manyfold: With larger source-samples, we will 
be able to derive a real chemical evolutionary sequence with one of the goals 
to use chemistry as an astrophysical clock. Furthermore, understanding the 
chemical differences is important to use the molecular fines as astrophysical 
tools to investigate the physical processes taking place. Moreover, another 
current hot topic is planet formation, and in this context astro-biology is a 
rising subject. In this regard understanding astro-chemistry and detecting 
new and more complex molecules in space is paving the way for future astro- 
biological science. 
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